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A STUDY OF Tri CHAIN-NELAXATION MATHOD OF NETWOKK ANALYSIS 


I. Introduction. 

(A) In the analysis of electrical circuits of the type commonly 
encountered in radio and communication devices, a problem of major 
importance is the choice of method by which the selution of a network 
is to be obtained. There are in general use today two principal methods, 
each having certain advanteges and disadvantages for the different types 
of network. A need has long existed for a single method which would be 
universally applicable, or nearly so, and which would have the added 
acvantage of reducing the effort required for solution. In a recent 
article (Bibliography 1), the "Chain-Relaxation" method has been pre- 
sented as a possible fulfillment of this need. It is the purpose of 
this paper to illustrate the application of this new method of analysis 
a8 a means of comparison with the more familiar procedures. 

(B) For the past half-century or more, there has existed e grow- 
ing interest in the mathematical analysis of electrical networks. As 
the field has grown wider and more complex, and circuits have become 
more complicated, the amount of calculation and effort needed for solu- 
tion of networks has increased many times. Today, there are in existe 
ence several electronic and electro~mechanical devices in laboratories 
and institutions wnich have been expressly designed for the purpose of 
analyzing the more involved circuits, Such machines can do the work 
of many individual mathesatieiens, who could conceivably spend hundreds 
Note: Because of time lost during the period when major effort would 
normally be applied to thesis work, resulting from a serious family 
ilineas, this study of the newly-introduced chain-relaxation method 


has been accepted as a master's thesis, in lieu of e more comprehensive 
project which had originally been chosen. 
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of man-hours on a single problem, in a matter of minutes. It is indeed 
fortunate that such advances have been made, and such analyzing mech- 
anisms will undoubtedly continue to be improved, but they are tremend- 
ously expensive items of equipment and require highly~specialized per- 
sonnel for their operation, Unless a problem is of rather great 
importance and wide applieation, it is not profitable to devise a 
computing mechanism for its solution and recourse must be made to clas- 
sical methods of computation. 

These classical methods of network analysis of lumped—perameter 
Circuits are usually designated as the "branch currents", "mesh-" or 
"loop~currents", and "nedal” methods. They have in common the technique 
of forming sets of simultaneous linear eauations on the basis of the 
application of certain slesentary laws of circuit behavior. These laws 
are: Kirchhoff's first rule, which states that the mm of the currents 
approaching and leaving a node mist equal zero, Kirchhoff's second rule, 
which states that the mum of the voltage drops around a closed path 
must be equal to zero, and Ohm's law, modified, which states that the 
current in any branch of a network is equal to the product of the ad- 
mittance of the branch and the net voltase across the branch. As 
normally written, these laws are: 

(a) 21 +0, at a node, 
(b) ZV = O, for a closed path, 
(e) ZV = I 
or 
21 = E, for a branch, 
where the definitions to be found in Guillemin, Bode, etc., (see Bibli- 


ography), for loop, node, and branch are intended. In these equations I 
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is the current in a branch, V is the potential of a node at which it 
terminates, and & is the electromotive force of any generator in the 
branch plus the electromotive forces produced in the branch by the 
action of vacuwa tubes or mutual inductances. 

The primary factor in the choice of method of attack of ea problem 
of other then extreme simplicity is the mumber of unlmowns that will 
have to be cetermined, Jor each of the three methods mentioned, cri- 
teria have been developed (libliography 2, 3, 4, 7) which permit de 
termination of the number of equations which can be expected to de- 
velop following the application of the method to a given network, If 
the number of branches is b, anc the number of nodes is n, then in the 
branch-currents method, there are n-] equations resulting from the ap-= 
plication of Kirchheff's first law to n-]) nodes plus b equations re~ 
sulting csrom the application of Chu's law, modified, for each branch, 
therefore b +n - 1 equations result, 

For the mesh-currents method, the resulting number of equations is 

- (mo and for the nodal method the number of equations is (n - 1). 
(A very good section en this toric is found in Bede's (Hibliography 3) 
first chapter.) 

It is evicent that the branch-current metnod will vive the greatest 
number of equations and it is for this reason that it has been almost 
entirely discarded except for very simple configurations, The mesh 
and nodal methods are those in general use today, anc ere the ones to 
be consicered here, 


_, the mesh-currents wethod was formulated by Haxwel] late in the 


last century and is perhaps more widely known than the nodal. method, 
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which has achieved popularity within comparatively recent years. 4ach 
has advantages end disadvantages. As an exmmple of the variation whieh 
can occur in the number of equations, Figure (1) of Plate I shows a 
cireult wicich by leop methods might require eight equations (assuming 
the branches shown could net be simplified) as against three by the 
nodal method, Figure (2) of the same plate sliows a circuit requiring 
two equations by the loop method as against seven by the nodal method. 
in general, it may be stated that where the nuxber of nodes is to re~ 
main fixed and branches added or taken away, ag in an experimental 
set-up, the nodal analysis is preferable, as the number of equations 
remains unchanged, requiring only a change in rows or colums of the 
eystom determinant if branches are added or taken away. Similarly, in 
& circuit in which portions of branches are tc be varied, keeping the 
number of meshes constant, the meai—current method is preferable. The 
nodal method is currently finding its widest use in analyzing vacuum 
tube circuits and in warp cases results in considerable reduction in 
the amount of labor required for solution. 

There does exist, however, one serious drewback to the nodal 
exelysis., The present  :rocedures cannot be applied if the network 
contains an arbitrarily located generator of negligible internal 
diapedance or having reguleted terminal voltage, as it becomes in- 
possible to replace the yvoltuge source oy an equivalent current 
source as is required by the nedal method. The manner of replace- 
ment is clearly presented in many textes (see, for example, Bibli- 
ography 3, 7). A second disadvantage lies in the fact that present 
sethods do not permit writing nodal equations directly from a circuit 


containing mitual inductance. Freliminary calculations are required. 
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A concise presentation of the manner of calculation is given in 
Gardner and Barnes, "Transients in Linear Systems", (Bibliography 4). 
It is worth mentioning at this point, that this second drawback is 
possibly on the way toward removal. Reed (Bibliegraphy 6) has worked 
out and presented a paper in which eae straightforward system of writing 
nodal equations directly frem the network, and also of taking into 
account the presence of voltage zenerate@s of negligible internal 
impedance, is outlined. if we assume that Reed's method is valid, 
then the mesh and nodal metheds may be considered of equal merit, and 
the choice of methods ia dictated by the network characteristics. As 
before, the method sheuld be chosen which involves the lesser amount 
of calewlation, or which more quickly obtains the desired reault. 

From the foregoing, it can be seen that having selected one of 
the two available methoda of analysis, and assuming that all possible 
use has been made of network simplification theorems prior to coumen~ 
cing the analysis, a certain fixed amount of computation will remain. 
if, therefore, another method of solution were developed which con- 
sistently gave fewer unknowns in the majority of problems of the type 
encountered today, its edontion should permit a freat reduction in the 
amount of effort reguired for solution of network problems. Such a 
method has been evolved and presented by Dr. L, Tasny-Tschiagsny, in 
the Journal of the Institute of Rlectrical Engineere (Liblicgraphy 1). 
li, Presentation of the Kethed. 

It is not the intention to here reproduce the article mentioned 
above, as it ig published and available, but an abstract of the salient 
features will be included in order that the illustrative exemples given 


te effect a partial comparison may have meaning. Let it suffice to say 
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that the author of the referenced article was inspired by Southwell 
(Bibliography 5), whe intredueed the method of "residuals" as the 
basis for his relaxation metiiods in stress calculations. The chain- 
relaxation method is an adaptation of these methods to use in electrical 
circuit work and provides a straightforward means of solution for most 
networks of the form commonly occurring in the field of radio and com 
munication as it is today. 

the chain-relaxatian method may be considered as a "forcing" 
method, in which as few as possible arbitrary branch currents are 
assumed and from them are derived node potentials, other branch cur- 
rents, driving currents of current generators, and e.m.f.'s of voltage 
generators by the wae of Kirchhoff's first rule and Ohn's lew, modified. 
Because of the arbitrary assumption of initial branch currents, it be~ 
Comes necessary at certain nodes to add "residual" currents in order te 
make the assumed current and potential distributions physically possible. 

A very few basie rules render the actual application of the method 
simple, if the meaber of generators is g ami the number of arbitrary 
aseumptions gs, then the number of residual currents (which are the un- 
knowns) will be given by: 

(d) ZBR=ZB-G- 

A second requirement of the method is that all branches contain- 
ing voltage sources be free of impedance, This can be accomplished 
in one of two ways ~ by arranging for two branches in series, one 
containing the senerator impedance and the other a cero impedance 
generator, or by replacing the voltage generator with its equivalent 
current generator having an impedance branch in parallel, ‘This is 


illustrated in Figure 3, Plate I. 
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Since it would be a remote chance that the asaimed currents were 
the actual values, the procedure is repeated s - 1 times with different 
arbitrary values and the results tabulated, These results will also be 
incorrect but by the principle of superposition the e.m.f.'s of genera 
tors or their driving currents may be made equal to their actual values 
and the residual currents may be made equal to zero, The values of the 
g multiplying factors will be unknown, but can be found by solving a 
set cf linear equations which are formed pronerly to fulfill the re- 
quirements. 

lil. The procedure is best illustrated by an exemple which will be 
taken step by sten, but, in general, one node is selected for reference 
and an arbitrary current assigned te a braneh terminating thereon. Uy 
successively applying Kirechhoff's first rule and Ohm's law, modified, 
adjacent node potentials and branch eurrents are derived until it be- 
comes necessary to wake an additional assumption before proceeding, 
Residual currents are assirned at any node whenever the currents neet- 
ing there do not comply with Kirchoff's first rule, 

Results are most corveniently set down in tabular form, Thera 
are four principal pointe in the procedure: 

(a) Schedule of steps. This is formed from the diagram of the 
circuit and indicates the direction in which the computation is to 
proceed, A symbol is umed to indicate steps of special neture, such 
as an arbitrary assumption, 

(b>) The cperations table. The schedule of steps is set down as 
Colum 1 of a table, in colum 2 is indicated by mathematical for- 
milae the nature of the computation for that ster, ami in s other 
colupms are placed the results of the computation for the g sets of 


arbitrary assumptions. 
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(c) The establishment of the equations. tere the results of the 
operation table are combined in such a way as to find the values of 
the s assumed currents that will render residual currents zero and 
e.m.f.'s or driving currents of generators equal to the know values. 
This step will be clearer when following an actual exemple, 

(d) The computation of values. The gs equations are solved for the 
necessary values of the originally assumed currents and from these val~ 
ues and the operations table, the values of the other branch currents 
and the node potentials are computed. The principle of superposition 
is employed. The expression for any quantity contained in a column 
of the operation's table is multiplied by tine correct value of the 
multiplying factor for that eolwm and the g products are added. 

Considerable simplification is achieved if the initially assumed 
currents are ail considered unity. That is, for the first assumption, 
consider the current. equal unity in thet colwm, and other currents 
zero, then in second coluwm, second current is one and others are 
zero, evc. This, however, is not a requisite, as any numericel value 
could be assigned, and the correct multiplying factor would be deter~ 
mined wien the equations were solved. 

As an illustrative example, the simple network of Flate Ii is 
selected and generalized impedances used in the interest of brevity. 
In practical applications, muwaserical values would be inserted at any 
point at which a saving im computation could be effected. Since the 
network has four independent nodes, requiring four uninowns, but may 
be solved with only three mesh currents, tie mesh analysis is chosen. 
The solution follows on succeeding plate for one of the mesh currents, 


wnich is in fact also a braneh current. 
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On plate (IV), the circuit is reproduced to permit relettering to 
show the sequence of steps by the chain-relaxation method, Table I of 
this plate is the operations table. Arbitrary assumptions are indica- 
ted by circles around the numeral, and steps leading to the establish- 
ment of the equations by squares. The computation will be carried 
through for one of the assumed currents which corresponds to that 
determined by mesh analysis, and it will be shown that the two are 
equal. Node 0 is taken as reference, The solution follows the plate 
showing circuit and the two currents are shown to be equal on follow- 
ing sheets, 

Step 1. Assign an arbitrary value of one te the current through 
element z¢. Insert 1 in colum 3, and zero in colum 4 of table. 

Step 2. The potential at 2 is the voltage rise above the refer- 
ence node. Insert 1 x Z¢ in colum 3 and 0 x L6 in colum 4. 

Step 3, No further movement can be made without a new assumption. 
Assume 1, » the current to node Q through Z,. 

Step 4. ‘The potential of node 4 is the voltage rise above refer- 
ence or 1,45. Insert OG x dn in colum 3 and 1 x Ky in column 4. 

The procedure is now apparent and will not be detailed. A few 
notes are in order, however, At step 6, we do not know that this 
relationship is true and might be inclined to add a residual current. 
sinee only one residual current will be required, r=s~g-1, it 
is more advisable to let these currents stand and wait to apply the 
residual at node 7, es in this way the potential of node 7 can be 
obtained from previously computed values simply. Had we added a 


residual at node 4, this simple reletionship would not hold since 
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the axpression for v4 would not contain only the assumed and derived 
values preceding but the reaidual as well, In step ll, we satisfy 
Kirehhof{'s first rule by the addition of the residual J,,, ané in 
step 12 we set up the finel condition required for a solution. The 
equations are then formed as shown on Plate (V). In order to obtain 
a solution the residual mast equal zero and V5.9 must equal known gen- 
erator voltage. The value of Iya in column 3 is multiplied by 4, ane 
added to the value for La in colum |, multiplied by L, and the sum 
equated to zero. Likewise the values of Vj, in columm 3 and J, are 
multiplied by I, and 1, respectively and their sum equated to the 
knewn E. The solution for iy follows and is later shown to be equal 
to the corresponding current obtained from the mesh analysis. The 
solution would be completed by solving for 1, and computing the values 
of any branch currents or node potentials from the values of coluwms 
3 and 4 as previously discussed by suming the respective products. 
I¥. Conclusions. 

it may appear that the saving is mall in this case. Another 
example is show on plates (YZ) and (21), and the initial equations 
are set up. This is e generalized version of ea twin-T network used 
by the author of the original article for illustrative purposes. Here 
it is apparent that either the mesh or nodal methods would require 
four unknowns, while the chain~relaxation method requires only two, 
and therefore is mach te be preferred. In general, it will be found 
that the number of unknowns is at least one Less than in either the 
mesh or nodal method when the circuit is of the type illustrated, 
The addition of vacuum tubes may but does not necessarily increase 


the number of unknowns in thig and the nodal method, while not 


ad 


Kervcrme hey Mersey wh) i thee ee Rie 1s sh rmderngee etd 
eben ow gt) eth eS Chee oe feelers il Be! elleeeye tendon 
ae yep See mY De A On Geis set) of Eiawete 
a keel & tel Sete cnt Adee haast o6f @ oe me UE ee 
ntedte of an ol (4) tee ot ce at ate) eat ree Rete 
eo ert Lame Sem 9? fn eeet Gees Spee Aether ey wheres ¥ 
ow A ye ahd we | emit nt ol te lee ape ey oe 
eae oe te Ss inde ele dee ley eet ud Delite 
te dee. cael 8 yl Se eek ort pales some of betas 
04 Of bodes om Les vas jlo some ps com oJ ro het inde 
CO 
wale ketlah wee Ot mt! se doe pblemyrsrten oSF of 
re Ae OM ee ee 
te 1k fey B wT otalAmy $0ed  bdierES aun Yu Se 
Te eee ot! a So ee eee a6 A See 


















tania 6 tT 


om we a 





Pees mart ee Coe gi eur ae 
: (23) |TV) my ee 

ee ome Behe 2 tit iter © One ope bee 
— men erase: at abt oD aes Ls eee ae ge 
erpet bee ee: toe & ae OS he ott bere 0b 
em ets ee Nee et eee alt eer emai“ et 
ee ee ee ee ee ee 
os ae oe ee ee ee ee Ge fe Geet ont Se 
DONTE me! et td tee oe ee ee wD 
cet hae fe kee ke eet ome, & ubie oeP 
oan OC ine peer tae 98) Me a A. oe Te seem om? 











affecting the mesh method, but this possible disadvantage is offset. 
by the fact that most vacuum tube circuits require more mesh than 
nodal equations. 

The addition of generators may work to the disadvantage of the 
chain-relaxation method, since the mwaber of unimowns can never be 
less than the number of generators. The other methods in certain 
networks require less unknowns as current sources or generators are 
added. 

For the types of circuits most frequently encountered, it appears 
that the chain-relaxetion method can be used to advantage. The con~ 
venient tabular form lends itself readily to computation and where it 
je desired to vary elements in the network, the effects of these varia~ 
tions can be readily determined, The form of the operation table per- 
mits the computation of only such values as are needed if the complete 
solution is not required, In any event, it should be possible to deter- 
mine by inspection of a circuit diagram prior to any computation whether 
er not the chain-relaxation procedure will reduce the required nuxber 
of unknowns below that of the better of the usual methods. 

It is this writer's opinion that the chain-relexation method is 
worthy of extensive trial ané further investigation. As yet, no formal 
set of rules has been made evailable covering such variables as the 
best sequence of steps to take or the most advantageous reference node 
to select. with practice, these points should present no great diffi- 
culty, as the optimum precedure will usually be apparent, but it would 
be of assistance to thome not familiar with the method. It is to be 
hoped that further work will be fortheoming in thia field which, at 
this time, holds at the very least the promise cf beceming an almost 


universally satisfactory method of network analysis. 
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